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Abstract 

O , Based on a sample of four million events collected by ALEPH from 1991 to 1995, 

a measurement of the forward-backward asymmetry in Z — > 66 decays using inclusive 
final states is presented. High-performance tagging of bb events in a wide angular 
^ . range is achieved using neural network techniques. An optimal hemisphere charge 

H ' estimator is built by merging primary and secondary vertex information, leading kaon 

identification and jet charge in a neural network. The average charge asymmetry, the 
flavour tagging efficiencies and mean 6-hemisphere charges are measured from data and 
used to extract the pole b asymmetry in the Standard Model 

A% = 0.1009 ± 0.0027 (siat) ± 0.0012 (syst) , 

corresponding to a value of the effective weak mixing angle of 

sin26'^ = 0.23193 ±0.00056 . 
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1 Introduction 

The measurements of the forward-backward asymmetry of 6-quarks from e'^e^ ^> Z ^ bb 
production offer the most precise determination of the weak mixing angle at LEP. This gives 
a sensitive test of Standard Model |jl| predictions of electroweak radiative corrections and 
hence constrain the allowed range of the Higgs boson mass. 

The polar angle distribution of the &-quark in the e~^e~ —>■ Z ^ bb process is forward- 
backward asymmetric: 

^'^bb 3 / 2/1 8 fc 

— (Tfeb - 1 + COS (Jb + -^FB COS &b 



d COS 9b 8 V 3 

where 6b is the angle between the incoming electron and the outgoing 6-quark. In this 
paper an improved measurement of Ap^ using inclusive 6-hadron decays is presented, based 
on the data sample recorded by ALEPH on the Z from 1991 to 1995. The analysis takes 
advantage of a reprocessing of LEP I data with improved charged particle tracking and of 
neural network techniques making maximum use of the available event information. 

The method of selecting events containing bb quark pairs is upgraded with respect to the 
previous analysis [Q. The 6-tagging algorithm based on the impact parameters of charged 
particle tracks is complemented with information from displaced secondary vertices, event 
shape variables and lepton identification. This leads to a 15% reduction of the statistical 
uncertainty. The estimate of the 6-quark direction is based upon the hemisphere charge 
method , but also incorporates information from fast kaon tagging and separated primary 
and secondary vertex charge estimators. The fraction of incorrect charge tags is reduced 
by 10% with respect to the previous measurement, implying a 20% further reduction of the 
statistical uncertainty. Another new feature is the control of systematic uncertainties by 
use of double tag methods for both flavour and charge tags, which yields a reduction of the 
systematic uncertainty by about a factor of two. 

Although lepton identification is used to select events in the 6-tag algorithm, no use is 
made in the hemisphere charge method of the lepton beyond that of an ordinary charged 
track in the detector. This reduces the correlation with the forthcoming ALEPH Apg analysis 
based on semileptonic final states. 

2 The method 

The measurement of Apg requires knowledge of the direction of the 6-quark from Z ^ bb 
decay. The quark-antiquark axis is estimated using the reconstructed thrust axis, the 
direction of which is given throughout this paper by a positive value for the cosine of the 
thrust polar angle, cos^. Each event is then divided into two hemispheres, F and B, by a 
plane perpendicular to the thrust axis. The forward hemisphere, F, is defined as the one 
into which the incoming electron points. The F-B orientation of the 6-quark is determined 
on a statistical basis by estimating the hemisphere charges. 

Hemisphere charges, Qp and Qb, are formed using a neural network designed to optimise 
the separation power between b- and 6-quarks. The neural network also performs well for 
lighter quark flavours, albeit not in an optimal way. A forward-backward asymmetry for 
flavour / at a given value of cos 6' is then proportional to the mean charge flow, (Qfb)' 



between forward and backward hemispheres in pure // events : 

{Qk) = {Qi-QL) 



-r {4{Qf - Qf) + 4(Q/ - Qf) 



where rip (rig) is the number of events with the primary quark from Z decay emitted into 
the forward (backward) hemisphere, and 6f = {Qf — Qf) is the average difference between 
the charges measured in the hemispheres of the quark and anti-quark, called the charge 
separation for flavour /. As shown in [Q, the same sample of events used to measure (Qfb) 
can also be used to extract Sf. This can be understood by considering a single hemisphere 
charge measurement, Qf, which can be written as : 



-j +nf and Qf- = -j +nf-, 



where TZ represents the measurement fluctuation due to fragmentation and detector effects. 
The product of the two hemisphere charges then averages to : 

{QfQf) = {QfQb) = ^ + (^/%) ' 

defining Sf = —Sf. The measurement fluctuation correlation, (TZfTlf), arises from effects 
of charge conservation, sharing a common event axis and crossover of particles close to the 
hemisphere boundary. It is then useful to define : 

5} = <y\QL) - <y\QL) 

= s) - ^iUfUf) - {Qi^f + {qL)^ 

= [Sf{l + kf)]\ (2) 

where Q(ot is the total charge measured in an // event. Thus, the observable Sf is equal to 
Sf to within a correction term, kf, depending on the polar angle and taking an average value 
of 9% for heavy flavours in this analysis. Compared with the charge separation itself, the kf 
correction term is less sensitive to the details of quark fragmentation and detector resolution 
0. Relying on the Monte Carlo prediction of kf, it is possible to extract Sf by fitting the 
quantities Sf to a range of measurements of the flavour combined S using the relation: 

^' = E ^/ (^/)' = -4(QfQb) - (Qfb)' + (Qtot)', (3) 

where Vf are flavour purities of the event sample under study. As described in Section ||, 
these purities are also measured in the data. 

For each value of cos 9 the mean charge flow, the purities and the charge separations are 
measured. The 6-quark asymmetry is then determined according to Equation (||), averaged 
over quark flavours: 

(QfbW) = E Vf{e)Sf{e)AUe). (4) 

f=u,d,s,c,b 



where 



Rearranging Equation (|) 



■^^"C)-^^!^- (^' 



AU<>) = jyskw ( **^™''''* " ^ 'PfiO)'>fW'^UO) j . (6) 

The 6-quark forward-backward asymmetry, Ap^, is extracted from a simultaneous fit to the 
measurements in nine angular bins in the range < cos0 < 0.95. A tiny, but non-zero, 
correction is applied to take into account the flavour dependence of the acceptance in each 
angular bin. 

3 The ALEPH detector 

The ALEPH detector is described in detail in P, and its performance in 0]. The tracking 
system consists of two layers of double-sided silicon vertex-detector (VDET), an inner 
tracking chamber (ITC) and a time projection chamber (TPC). The VDET single hit 
resolution is 12 /xm at normal incidence for both the r0 and rz projections and 22 /im 
at maximum polar angle. The polar angle coverage of the inner and outer layers are 
I cos 6*1 < 0.84 and | cos 6*1 < 0.69 respectively. The ITC provides up to 8 r0 hits at radii 
16 to 26 cm relative to the beam with an average resolution of 150 /im and has an angular 
coverage of | cos 9\ < 0.97. The TPC measures up to 21 space points per track at radii 
between 38 and 171cm, with an rcj) resolution of 170 /im and a z resolution of 740 /im and 
with an angular coverage of | cos 6\ < 0.96. In addition, the TPC wire planes provide up to 
338 samples of ionisation energy loss (dE/dx). 

Tracks are reconstructed using the TPC, ITC and VDET which are immersed in a 1.5T 
axial magnetic field. This provides a transverse momentum resolution of a{l/pT) = 0.0006 
(GeV/c)~^ for 45 GeV muons. Multiple scattering dominates at low momentum and adds a 
constant term of 0.005 to a{pT)/pT- 

Outside of the TPC, the electromagnetic calorimeter (ECAL) consists of 45 layers of 
lead interleaved with proportional wire chambers. The ECAL is used to identify photons 



and electrons and gives an energy resolution a{E)/E = 0.18/JE/GeY + 0.009 for isolated 
particles. The hadron calorimeter (HCAL) is formed by the iron of the magnet return yoke 
interleaved with 23 layers of streamer tubes. It is used to measure hadronic energy and, 
together with two surrounding layers of muon chambers, to identify muons. 

The information from the subdetectors is combined in an energy flow algorithm ^ which 
gives a list of charged and neutral track momenta. 

Recently the LEP I data have been reprocessed using improved reconstruction algorithms. 
In particular, the VDET hits are distributed among tracks according to a global x^ 
minimisation procedure which improves the hit association efficiency by more than 2%. 
Information from TPC wires is now used in addition to the pad information to improve the 
coordinate resolution by a factor of two in z, and by 30% in r0 for low momentum tracks. 
Similarly the pad information is added to the dE/dx information from TPC wires, increasing 
the fraction of tracks with a useful dE/dx measurement to almost 100% and providing on 
average a two sigma separation between pions and kaons with momenta above 2.5 GeV/c. 



4 Monte Carlo simulation 



The analysis makes use of a Monte Carlo sample (MC) of 8.1 million simulated liadronic Z 
decays as well as two dedicated heavy flavour samples of 4.9 million Z ^ bb decays and 2.4 
million Z —>■ cc decays. The simulation is based on JETSET p| with parameters tuned to 
reproduce inclusive particle spectra and event shape distributions measured in hadronic Z 
decays ^. 

In the MC simulation, the most relevant physics input parameters have been adjusted, 
using the re- weighting technique, to the recently measured values |l^, |ll| shown in Table 



The fragmentation of heavy quarks into hadrons is assumed to follow the model of Peterson 



et al. [0 with parameters tuned to match the values of the mean heavy hadron fractions of 
the beam energy reported in the table. 



Physics parameter 


World average value 


< Xb > beam energy fraction 


0.702 ±0.008 


< Xc> beam energy fraction 


0.484 ±0.008 


rich in &-hadron decay (i^° and A inch) 


5.44 ±0.09 


Bg fraction 


0.100 ±0.012 


Afc fraction 


0.099 ±0.017 


B+ hfetime 


1.656 ± 0.025 ps 


B^ lifetime 


1.562 ± 0.029 ps 


Bg lifetime 


1.464 ± 0.057 ps 


6-baryon lifetime 


1.208 ± 0.051 ps 


g ^ bb rate 


0.00251 ± 0.00063 


g ^ cc rate 


0.0319 ±0.0046 



Table 1: List of physics input parameters [|T0|, |TT| to the Monte Carlo simulation which are 
used for re-weighting. 



5 The neural net 6-tag 

Due to the long lifetime and high mass of 6-hadrons, 6-jets have several characteristic 
properties. Six discriminating variables are combined using a neural network to tag 6-quark 
jets. A similar scheme is used by ALEPH to identify 6-quark jets in searches for neutral 
Higgs bosons conducted at LEP II [p^ . 

The jets are clustered with the JADE algorithm using a ycut value of 0.02, and for each 
jet six variables are defined: two of them are lifetime-based; a third one is based on the 
transverse momentum of identified leptons; the last three are based on jet-shape properties. 
These quantities are 

1. Pjet: probability of the jet being a light quark (uds) jet based upon impact parameters 



of tracks in the jet |14 



2. Axsvy^'- the vertex x^ difference between assigning tracks in the jet both to the secondary 
and primary vertices compared to assigning all tracks to the primary vertex. This 
is based upon a secondary vertex pattern recognition algorithm which searches for 
displaced vertices via a three-dimensional grid point search |]T3[ ; 



xlO 



4000 



ALEPH 



o Data 
- MC 



1.075 

1.05 

1.025 

1 

0.975 

0.95 

0.925 

0.9 



V 



Data/MC 



+, 




V.,_,„..y^V«V*• Vi 



\i 111! TJilTTii ' 



W ^ HizYt ^' ^H 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

b-tag 



Figure 1: Distribution of the neural net hemisphere 6-tags and the ratio between data and 
MC. The shaded region is the distribution for 6- hemispheres in the MC. 



3. pt of identified leptons [0 with respect to the jet axis; 

4. S^: the boosted sphericity of the jet defined to be the sphericity of final state charged 
and neutral particles in the rest frame of the jet; 

5. Multiplicity/ln(Ejet): the charged and neutral particle multiplicity of the jet divided 
by the natural logarithm of the jet energy in GeV; 

6. Sp|.: the sum of squared transverse momenta, p^, of each charged or neutral particle 
with respect to the jet axis. 

Among the six input variables, the lifetime tags have the largest weight. However, the 
inclusion of the other variables increases the 6-quark discriminating power, especially at the 
most forward angles that are not covered by the vertex detector. A hemisphere 6-tag is 
defined as the maximum neural net output value among the jets in the hemisphere with 
energies exceeding 10 GeV. 



Prior to the application of the hfetime reconstruction algorithms in the MC, the track 
parameters have been given additional smearing and the VDET additional hit inefficiency 



according to a procedure described in |^. This is necessary in order to render the MC 
distributions of the hfetime tags in good agreement with data. Smah discrepancies still 
remain and combine in the neural network to give the differences between the distributions 
of the 6-tag shown in Figure |l[ Deviations up to about 5% between data and simulation 
are seen in some 6-tag bins and for this reason the 6-quark purity as a function of 6-tag is 
determined from the data themselves. 

6 The neural net charge tag 

A neural network is also used to determine hemisphere charges. This technique has been 
employed earlier for CP violation studies in B^/B^ decays as described in [^ . In the present 



analysis information regarding lepton identification is left out from the tag. 

In order to achieve high tagging performance independently of the 6-hadron species, 
momentum and polar angle, several charge estimators are combined. These charge tags are 
complemented with other variables, which do not carry information about the 6-hadron 
charge in themselves, but are correlated with the relative tagging power of the charge 
estimators, thereby allowing an optimal combination to be achieved in the whole phase 



space. The two types of input variables are described separately in Sections |6]T| and |0| . 



6.1 Charge estimators 

Eight charge estimators are used to provide information on the initial charge state of 6-quark 
hemispheres as inclusively as possible. These are: 

1. The jet charge. This is the charge estimator used in the previous analysis 0: the 
weighted sum of particle charges in a hemisphere, where the weights are the particle 
momenta along the thrust axis, pl, raised to the power k. Four different values 
of K are used for inputs to the neural net: 0.0, 0.5, 1.1 and 2.2. These focus on 
correspondingly higher values of track momenta but are of course highly correlated. 
Such charge estimators work well for all species of 6-hadrons. 

2. The secondary vertex charge. Using a topological vertexing algorithm combining 



information from all charged tracks in the hemisphere |T8[, a best estimate for a 
secondary vertex position is obtained in each hemisphere. An estimate of the charge 
of this vertex is calculated as: 

Qvtx = "^ WiQi, (7) 

i=tracks 



where Wi, calculated by a dedicated neural net ||T8[, is the probability that track i comes 
from the secondary vertex and qi is the track charge. This estimator provides a high 
quality charge tag for B^ hemispheres, and also helps to indicate which hemispheres 
are more likely to contain a neutral B meson. 

3. The weighted primary vertex charge. For B^ hemispheres the secondary vertex charge 
carries little information, but in this case the fragmentation tracks close in phase space 



to the B^ have some correlation with the primary 6-quark charge. Therefore a charge 
estimator, Qp^tx, is calculated according to : 

Qpvtx = Y. (l-^i)PL9i/ E (i-Wi)pL- (8) 

i=tracks i=tracks 

The weighted secondary vertex charge, Qsvtx, is calculated in a similar manner, 
replacing the weight (1 — Wi) pl^ by Wi (pl) with the aim of improving the tagging 
for both charged and neutral 5-decays via leading particle effects from the secondary 
decay. 



5. Fast kaon identification is formed by another dedicated sub-net [|18[ trained to identify 
charged kaons from 6-hadron decays. The output for the most kaon-like particle in the 
hemisphere is signed by the charge of the particle and used as a charge estimator. 

6.2 Control variables 

The following variables are used as inputs to the neural net in order to provide it with some 
topological and 6-hadron specific information: 

1 . I cos 6 1 is included since it is correlated with the quality of the secondary vertex 
reconstruction and hence with the tagging power of the estimators relying on the 
separation between tracks from the primary and secondary vertex. 

2. The reconstructed b-hadron momentum is included because the relative accuracy of 
the various charge estimators depends on the 6-hadron momentum. An estimator of 
this momentum is constructed from the jet closest to the line-of-flight of the 6-hadron. 
The estimator is the sum of the charged track momenta in the jet, weighted by the 
probability that they come from the secondary vertex, and the projections of the 
neutral momenta onto the line-of-flight of the 6-hadron. The missing energy in the 
hemisphere is also added. 

3. The reconstructed proper time of the b-hadron is used based on the reconstructed 
6-hadron momentum and the measured decay length. The intention here is to 
incorporate the increased probability of B^ mixing at long proper times. 

4. The spread of track separation weights i.e. the width of the Wi factors distribution in 
a given hemisphere. This allows the net to de-weight those charge estimators which 
suffer from an ambiguous allocation of tracks to the primary and secondary vertices in 
cases of high charged multiplicities and/or poor vertexing. 

For the asymmetry calculation, the difference between the neural net outputs of the two 
hemispheres is used. This is shown in Figure ^ together with the MC prediction. The charge 
separation is seen to be reasonably well simulated, but with a slightly larger width in data 
than in MC (there is also a small difference in the average value of the two distributions, 
but this is too small to be visible on the plot). 



16000 


- 








14000 
12000 


• Data 
— MC 


/^ 


\ 


ALEPH 


10000 
8000 


— 


'^^^'o 


X 




6000 


r jy 




^^o 


\ 


4000 


r jT 




. 


V 


2000 
n 











-2 -1.5 -1 -0.5 0.5 1 1.5 2 

Qfb 



Data/MC 




Figure 2: The charge difference between forward and backward hemispheres measured in 
the selected event sample labelled E in Figure ^j. The MC distribution corresponds to an 
asymmetry of Ap^ = 0.0967. 



7 Event selection 



The data set used for this analysis consists of four million hadronic Z decays recorded by 
ALEPH during the period 1991 to 1995 in a centre-of-mass energy range of Mz ± 3 GeV. 
Events are selected according to the standard ALEPH hadronic event selection based on 
charged track information |T^ . This selection has an efficiency of 97.5 % and the backgrounds 
from Z decays to t^t~ and 77 interactions are estimated to be ~ 0.3% each. These 
backgrounds are reduced to the 10~^ level after the application of the subsequent cuts and 
can be safely neglected. 

The average beam-spot position is determined every 75 events and used to constrain the 
event-by-event interaction point. In order to remove lifetime information from the tracks, 
these are projected onto the plane perpendicular to their parent jets (unless they point 
behind the interaction point found in a ffist iteration). Combining these projections with 
the beam-spot position fixes the interaction point to a precision of 50 x 10 x 60 /xm^ in 
horizontal, vertical and beam directions respectively. 
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Figure 3: Upper plots: Subdivision into 6-tag intervals with numbers above the plots 
enumerating the intervals. The shaded histogram is the expected 6-tag distribution for 
6-hemispheres, shown separately for two regions of polar angle. Lower plots: The event 
acceptance, shown shaded and labelled E, in the plane of the two hemisphere 6-tags. 



The thrust axis is determined in each event from all charged tracks and calorimeter 
clusters using the ALEPH energy flow algorithm [|^. To ensure that a good fraction of the 
event is contained within the detector acceptance, the cosine of the thrust axis polar angle 
is required to be less than 0.95. In order to be able to evaluate the hemisphere 6-tag and 
charge tag it is required that each hemisphere has at least one charged track and at least 
one jet with energy greater than 10 GeV. With these selection cuts a sample of 3,734,425 
hadronic Z decays is retained for the remainder of the analysis. 



In order to make use of double-tag methods for determining flavour purities, the 6-tag 
range between zero and one is divided into eight sub-intervals as shown in Figure 0. The 
672,111 events which are located in the region indicated by the black area in Figure]^, labelled 
E, are the ones accepted for the ApB analysis. The acceptance is less strict in the forward 
region which is not covered by the vertex detector and therefore has a broader distribution 
of 6-tags. 

8 Determination of flavour purities 

The fractions of primary quark flavours in the selected event sample are measured in the 
data using the double tag method, that has been used previously in e.g. the ALEPH Rj, 
measurement [0, |2^ . 

In each of nine angular intervals, eight 6-tag sub-intervals are defined as indicated in 
Figure 0, out of which only seven are mutually independent and the eighth is considered as 
"not tagged" . These give rise to 28 doubly tagged event classes and 7 singly tagged classes: 
riij and riis, which are related to the tagging efficiencies by: 



Hij 



f 

n,s = E2efi?;(l-5]e|4.)^had; ^<7, (9) 

/ i<7 

where e{ is the probability for a primary quark of flavour / to be located in interval i, c(, 
describes correlations between the 6-tag values in the two hemispheres, Rf is the Z ^ ff 
fraction of the Z hadronic decays and rihad is the total number of events in the sample under 
consideration. The flavour index, /, runs over three flavour classes only [uds, c and b), hence 
efficiencies and correlation factors are averaged over the three light flavours. 

Various assumptions are made to reduce the number of unknown parameters in 
Equation @. The hemisphere correlations are taken from MC simulation and the values of 
Rf are taken from previous measurements. Furthermore, it is necessary to flx the ratio of the 
two smallest efficiencies, where e"'^^ ~ 0.0008 and £3'^'' ~ 0.0013 according to the simulation. 
The ratio between these numbers is taken from simulation. 

The remaining 20 unknown efficiencies ef are determined from a fit using Equations @ 
leaving 15 remaining degrees of freedom in the flt. 

The fltted efficiencies are combined to give event tag efficiencies by summing over the 
fiducial region labelled E in Figure ^ 



/ - sr J J J 



Ee/ej4, (10) 



and from these efficiencies the three purities are constructed as: 

-pf, = ^_ERf_ ^ ^^^) 

e^Rb + e%Rc + e'%'^{l — Rb — Re) 

where the used Rf fractions are shown in Table ^. The result of the fit is shown in Table ^ 
The covariance matrix of the fitted efficiencies is propagated to the event-tag purities with 
diagonal elements also shown in Table ^ and correlation coefficients being typically —90% , 
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Fitted purities 




Fit-MC 


cos(^) 


uds 


c 


h 


uds 


c 


h 


0.0-0.1 


157 ±46 ±3 


1237 ±69 ±7 


8606 ± 34 ± 15 


-33 


96 


-63 


0.1-0.2 


239 ± 56 ± 5 


1170 ±77 ±7 


8591 ±32 ±14 


43 


6 


-49 


0.2-0.3 


208 ± 48 ± 4 


1156 ±69 ±7 


8636 ±31 ±14 


15 


13 


-28 


0.3-0.4 


246 ± 55 ± 5 


1098 ± 74 ± 6 


8657 ±31 ±14 


64 


-26 


-38 


0.4-0.5 


200 ± 36 d= 4 


1064 ±56 ±6 


8736 ± 40 ± 14 


32 


8 


-40 


0.5-0.6 


147 ± 42 ± 3 


992 ± 63 ± 6 


8863 ±31 ±14 


-15 


22 


-7 


0.6-0.7 


197 ±52 ±5 


809 ± 73 ± 5 


8994 ± 29 ± 14 


42 


14 


-56 


0.7-0.8 


422 ± 94 ± 16 


366 ±109 ±3 


9212 ±30 ±14 


238 


-242 


-4 


0.8 -.95 


405 ± 203 ± 5 


883 ± 250 ± 8 


8721 ±57 ±14 


27 


342 


-369 



Table 2: Flavour purities of the event sample labelled E in Figure |^ as obtained from the 
double tag fit. The values are given in units of 10~^. The second error is from MC statistics 
contributing through the assumed hemisphere tag correlations. 



±45% and —80% for the uds-c, the uds-b and the c-b purity correlations, respectively. The 
X^ of the nine independent fits at different cos 6 averages 0.99 per degree of freedom. 

Already from Figure |l| it is seen that the efficiencies in the highest 6-tag bins are higher 
in the data than in the MC. The fit determines that this enhancement is slightly larger for 
light and charm quarks than for 6-quarks, resulting in a 6 purity which is about 0.5% lower 
than predicted by the MC in the selected sample. 

9 Determination of quark charge separations 

The mean values of the hemisphere charge separations, Sf, are determined from data in each 
bin of cos 6 using the following procedure. 

The data sample is subdivided into 14 bins with fiavour compositions ranging from almost 
pure uds to almost pure b fiavour. The sum of squares of the two hemisphere 6-tags is used 
to define the bins and the fiavour compositions of these new bins are derived from the fit to 
the data described in the previous section. 

In each bin of this event 6-tag variable, 6 is measured according to Equation @. A fit is 
performed using the 14 6 measurements and the right side of Equation (^ with the fiavour 
specific values 6b, Sc and 6uds as free parameters. Since the 6-tags in the 14 bins bias the 
charge estimator differently, additional assumptions taken from MC are needed in order to 
compare the 14 6 measurements. For each fiavour, the differences between the 5 in a given 
bin and the 6 in the E-region, jS = 6e — S^n, are held fixed at their MC values. Thus, for 
each fiavour, the fit shifts 5 by a fixed amount in all bins, minimising the following x^'- 



X 



14 

E 
j=i 



6^ 

"meas,j 



E%^Vi{6frf-(3l 



I ^^^incas.i 



:i2) 



where i labels the 6-tag bins and / the fiavours (6, c and uds). The 6f appearing in the x^ 
refers to MC events in the E-region and the correction factors Vf are the three free parameters 
of the fit. 
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The inputs to the fit, averaged over all cos^, are shown in Table |^ for b- and c-quarks 
together with the flavour combined output and the MC prediction. The fit reproduces 
accurately the 6-tag dependence of 6. The angular dependence of the fitted parameters is 
shown in Figure || and, since the ratios between fitted and predicted values are consistent 
with being constant over angles, their average values are quoted below: 

rb = 0.995 ±0.004, 

rc = 1.031 ±0.011, 

ruds = 1.009 ± 0.009 . 



Vc 


n 


Pc 


Pb 


Observed 
6 


Expected 
6 


Fitted 
6 


5.2 


0.3 


0.0676 


0.1650 


0.4413 ±0.0061 


0.4407 


0.4430 


8.6 


0.8 


0.0808 


0.1559 


0.3877 ±0.0039 


0.3851 


0.3878 


11.8 


1.6 


0.1016 


0.1663 


0.3347 ±0.0036 


0.3341 


0.3371 


15.5 


3.0 


0.1119 


0.1493 


0.3115 ±0.0035 


0.3083 


0.3124 


20.0 


5.0 


0.0979 


0.1496 


0.3039 ± 0.0033 


0.3003 


0.3042 


24.9 


8.2 


0.0884 


0.1284 


0.3030 ±0.0031 


0.2960 


0.3001 


29.7 


12.2 


0.0729 


0.1186 


0.3048 ± 0.0030 


0.2982 


0.3029 


34.4 


18.4 


0.0572 


0.1118 


0.3011 ±0.0028 


0.2994 


0.3026 


38.4 


27.7 


0.0465 


0.0945 


0.3080 ± 0.0026 


0.3028 


0.3072 


36.2 


45.3 


0.0237 


0.0813 


0.3144 ±0.0023 


0.3092 


0.3130 


17.4 


77.7 


0.0098 


0.0445 


0.3113 ±0.0021 


0.3146 


0.3157 


12.5 


85.8 


-0.0104 


0.0290 


0.3276 ± 0.0028 


0.3272 


0.3275 


7.0 


92.5 


-0.0252 


0.0091 


0.3411 ±0.0029 


0.3421 


0.3420 


1.2 


98.8 


-0.0802 


-0.0455 


0.3910 ±0.0019 


0.3941 


0.3886 



Table 3: From left to right: bin-by-bin flavour purities (V in percent), biases (/? = 6e — <^bin), 
measured, MC expected and fitted values of 6un- The numbers in this table are integrated 
over all polar angles. 

The average correlation coefficients of the error matrix are —64% between uds and c, 
25% between uds and b and —49% between c and b. The average x^ probability is 0.72 over 
the independent fits in nine bins of polar angle. 

The fitted values of 6f are translated into the charge separations, 6f by Equation (H), 
where the charge correlation correction factors, (1 ± kj), are taken from MC. These factors 
vary with polar angle, averaging 1.086 for c-quarks and 1.092 for 6-quarks. 

In the fit of Ap^ from Equation (|^), the values for the charge separations, Sf{9), are thus 
the MC values corrected by the rf{6) factors. The ratios between individual light quark 
charge separations are taken from MC. In summary, the following values are used (averaging 
over polar angle and showing only the diagonal elements of the statistical error matrix): 



St 



-0.3143 ±0.0014, 
±0.3733 ±0.0039, 
-0.1929 ±0.0016, 
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Figure 4: The fitted values of 6b and 6c as a function of polar angle together with the values 
predicted by the MC. 

6a = -0.1547 ±0.0013, 
6u = +0.2598 ± 0.0022 . 



This charge separation corresponds to an average event mis-tag of 25.6% for bb events. The 
result of the asymmetry fit to the peak data is 

A^3(^ = 91.232 GeV) = 0.0997 ± 0.0027(stat). 

10 Systematic errors 

The determinations of the sample purity, the charge separation and the forward-backward 
asymmetry itself all rely to some extent on samples of Monte Carlo events. The features of 
the simulation most relevant for the analysis have been identified and varied within their 
errors in order to estimate the systematic uncertainty of the measurement as described in 
detail in the following sections. 
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10.1 The uds purity at large 6-tags 

In the fit yielding the flavour purities of the selected event sample it was assumed that ratio 
of the tiny efficiencies for selecting a light flavour with 6-tag > 0.98 and for selecting a light 
flavour with 0.94 < fe-tag < 0.98 was correctly predicted by the MC to be 0.62, averaged 
over polar angle. A systematic uncertainty on this ratio of ±30% is estimated by varying 
the track parameter smearing and the gluon splitting rates to heavy flavours in the MC. 
Propagating this uncertainty to the asymmetry yields a variation of ±0.00037. 

10.2 Correlations between hemisphere 6-tags 

The correlation factors (q- in Equation (^) are taken from Monte Carlo simulation. The 
influence of the correlations on the purities, via the fitted flavour efficiencies, is limited by 
Equation (0) because the 6-quark purity is high, 88% on the average. The correlation 
coefficient with the largest impact on the b purity is the negative correlation for observing 
both hemispheres in the highest 6-tag bin: c\i — 1 = —0.02 in the central region of the 
detector. It is due to the pull on the shared reconstructed primary vertex exerted by a 
secondary vertex from 6-hadron decay. If that coefficient is changed by an amount ±0.01, 
the b purity changes by the amount ±2.3 ■ lO"'* in the central region. The other coefficients 
carry impacts on the flavour purities of size lO"'^ or less for a one percent change in the 
coefficient. 

In order to estimate the uncertainty in the correlation coefficients, the input parameters 
to the MC given in Table are changed, one by one by the re-weighting technique, to the 
one a higher values. Two additional re-weighted MC samples are studied in order to cover 
effects of QCD gluon radiation and of fragmentation. In one sample the angles separating 
the two selected 6-jets are forced to agree with data and in another sample the charged 
multiplicity of the primary vertex is forced to agree with data. The largest effects on the b 
purities arise from changing the Xf, beam energy fraction, the charge particle multiplicities 
and the inter-jet angles, ranging from typically ±10~^ in the central region to ±10"'^ at the 
most forward angles. The individual changes are added in quadrature to produce the total 
variation in the asymmetry of 0.00001. 

There could be effects other than 6-physics parameters causing an inaccurate simulation 
of the correlation factors. The possibilities include inaccuracies in the simulation of vertexing 
and of QCD effects. Therefore, a study of additional systematic errors is carried out along 



the lines of the ALEPH Rh measurement |T^, [2^]. Here, the contribution to the correlation 
from a given variable v is calculated as: 

''^^ 2 iJ Mv)dv) iJ f,{v)dv) ^''^ 

where fi{v) is the fraction of hemispheres, among all hemispheres with the variable situated 
between v and v + dv, having the tag i. Similarly gi{v) is the corresponding fraction with 
the opposite hemisphere tagged by tag i. The chosen variables, v, are the thrust polar and 
azimuthal angles, an estimator of the 6-hadron momentum and the six input variables to 
the 6-tagging neural network. Since the most important correlations are those involving bb 
events, the ten most signiflcant of those are studied. In order to avoid having to subtract 
a very large non-6 background in some bins, a soft 6-tag is placed before evaluating the 
integrals. This soft 6-tag leaves a 64% pure 6-sample. 
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For the azimuthal angle, where no significant correlation is expected, the result of 
Equation (|13|) is indeed unity within a precision of 10~^. The rest of the variables do 
produce correlations, but the differences between data and MC are all between 10~^ and 
10~^ and hence their contribution to the error on the measured asymmetry (via the flavour 
purities) is again very small, in total 0.00001. 

Until now only correlations in bb events have been studied. Although the individual 
correlations for the lighter quarks have very little impact, the possibility of a combined 
effect of a large number of correlations has been investigated by simply ignoring all of them 
in the fit (except at the very most forward angle where geometrical effects are important). 
This results in a change in Ap-^ of +0.00002. A summary of all errors connected with the 
estimation of flavour purities is given in Table HI. 



Statistical error in the purity determination 
Light quark 6-tag efficiencies 
Hemisphere 6-tag correlations 



Total systematic error due to purities 



±0.00039 
±0.00037 
±0.00003 



±0.00054 



Table 4: Summary of systematic errors on the Apg measurement from the determination 
of flavour purities. 



10.3 Systematic uncertainties on charge separations 

Various assumptions contribute to systematic uncertainties on the determination of the 
charge separations, 6f. In a first step, Sb, Sc and Suds are extracted from a fit to data, where 
the assumed flavour purities and charge bias corrections in the different 6-tag bins contribute 
to the systematic error together with the statistical errors of the fit. In a second step, the 
constant of proportionality, (l + kf), between Sf and 6f contributes to the uncertainty. These 
systematic uncertainties are summarised in Table ^ and discussed in detail in the following. 

• Flavour purities - This is evaluated by using the purities predicted by the MC instead 
of the fitted purities. The full changes in the fitted values of Sf are propagated to the 
error on the asymmetry. 



Corrections for b-tag bias 
uncertainty: 



Two overlapping procedures are used to estimate this 



1. Each of the 6-physics parameters listed in Table |l| are increased one at the time by 
one standard deviation, using re-weighted versions of the MC, and the differences 
in the fitted Sb and Sc are propagated to the asymmetry and entered as systematic 
errors in Table |^ . These estimates include the effect of the parameters on the 
charge correlation corrections discussed below. As in the case of hemisphere 
tag correlations, the effect of re-weighting the MC to agree with the observed 
primary vertex multiplicity and the observed distribution of inter-jet angles are 
also included. Finally, the Vector /{Vector ± Pseudo scalar) ratio in D-meson 
production is varied within its experimental error of 0.05 [pl[]. 
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2. The values of Sc and Suds are fixed at their MC values instead of being left free 
in the fit. In this case Equation (^) can be solved directly for 61, in the E-region 
without any assumptions concerning bias corrections. The full difference between 
the result of this and the default result is taken as an independent contribution 
to the systematic error. This corresponds to an excursion in the fitted value of 6c 
by almost two times the other uncertainties on this quantity. 

Charge correlation corrections - One contribution to the systematic error from the 
(1 + kf) factors is found by varying the input parameters of the MC as described 
above. 

Another contribution is estimated, like in the case of hemisphere 6-tag correlations, 
from the difference between data and MC when building on the E-sample a "projection" 
of the hemisphere charge correlations along a series of observables, as follows: 

"■'*•' = m^)m ■ <"' 

where Q^{v) is the average positive or negative hemisphere charge, given the variable 
V in either the same or the opposite hemisphere, and P{v) is the probability density of 
V. Eight variables are chosen that are weakly correlated while covering the variables 
relevant for the charge determination. These are the thrust, the thrust axis direction, 
some inputs to the charge tagging neural net, the 5-tag and the 6-hadron momentum. 
The differences in (1 + k^) between data and MC are propagated to Table ^. 

This procedure relates to the simulation of charge correlations for events in the selected 
sample consisting mostly of bb events. Previous studies using the jet-charge technique 
[|, 1^ have assessed the uncertainty in kc separately. Using a weighting power of 
K = 0.5 the value kc = 0.085 ± 0.022 was found, where the error includes systematic 
uncertainties. This value is equal to the charm charge correlation in the present study 
and its error is propagated as an additional contribution to Table ^. 

Light quark charge separations - The light quark charge separations, 6u, Sd and 6s, are 
taken from Monte Carlo, but corrected by the overall correction factor from the fit. 
Systematic uncertainties are ascribed according to the results of which represent 
relative uncertainties of 2.1%, 4.0% and 1.5% for u, d and s respectively. These errors 
are much larger than the error on the combined light quark 5uds from the fit. However, 
when propagated through to A^^-q, the systematic uncertainties arising from light quark 
charges are small, in total 0.00003. 

Track parameter smearing - The extra smearing and hit inefficiency mentioned in 
Section ^ also affects the charge tag. The effect of dropping these corrections to the 
MC altogether results in a change in the asymmetry of +0.00016 which is taken as an 
additional contribution to the systematic error. 



10.4 Experimental systematic errors on Qfb 

For the measuring apparatus itself to produce a fake Qybi a forward-backward asymmetric 
bias in the charge measurement is needed. Such an effect could come from an asymmetry in 

16 



Source of systematic uncertainty 



Statistics in the 6f determination 



flavour purities 



< Xb > beam energy fraction 

< Xc> beam energy fraction 
rich in 6-hadron decay 

Bs fraction 
Ab fraction 
B~^ lifetime 
B^ lifetime 
Bg lifetime 
fe-baryon lifetime 
g ^ bb rate 
g ^ cc rate 
primary multiplicity 
inter-jet angles 
for charm 



v+p 



fixed Sc and S, 



uds 



charge correlations due to thrust 

charge correlations due to cos^ 

charge correlations due to 

charge correlations due to 6-hadron momentum 

charge correlations due to 6-tag 

charge correlations due to jet-charge 

charge correlations due to primary vertex charge 

charge correlations due to secondary vertex charge 

charge correlations due to kaon charge 

charm quark charge correlations 



light quark charge separations 



track parameter smearing 



Total systematic error due to charge separations 



^AU 



0.00074 



0.00015 



0.00007 
0.00000 
0.00003 
0.00004 
0.00004 
0.00003 
0.00002 
0.00002 
0.00003 
0.00008 
0.00004 
0.00001 
0.00004 
0.00019 



0.00037 



0.00001 
0.00008 
0.00007 
0.00002 
0.00000 
0.00012 
0.00010 
0.00011 
0.00000 
0.00016 



0.00003 



0.00016 



0.00093 



Table 5: Summary of systematic error contributions to the Apg measurement from the 
determination of charge separations. 



the detector material, since nuclear interactions of the produced particles with this material 
give rise to an excess of charge which is measured by the total charge, Qtot, in data 0. The 
asymmetry of the material is estimated by measuring the forward-backward asymmetry of 
photons that have converted to e^e~ pairs in the detector. This asymmetry is multiplied by 
the material charge component at each angle and subtracted from Qfb, causing a shift in 
ApB by 0.00016 ± 0.00032 which is taken as a contribution to the systematic error. 

Another possible source of error is the modelling of the magnetic field in the extreme 
forward region. This is known to cause a slight difference in the momenta measured for 
positive and negative particles at momenta close to the beam energy and a correction, which 
is not quite forward-backward symmetric, is applied as default on the data to take out this 
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effect. However, the effect on Apg from dropping the correction is very small (0.00002) and 
this is added to the error. 

10.5 Electroweak observables 

Assumptions regarding electroweak observables other than Ap-^, e.g. Rb and ^fb; belong to 
a special class, since they are important for the Standard Model interpretation of the result. 
Their values used in this analysis are taken from the Standard Model and are listed in Table ^ 
together with their impacts on the measured value of Ap-g. In case of the branching fractions 
Rf, the variation in one flavour is compensated by changes in the light quark R/s (/ = uds) 
according to their relative magnitudes. 



Parameter P 


v^ = 89.416 GeV 


v/i = 91.232 GeV 


v/i = 92.945 GeV 


M^b/5 P 


Ru 




0.17211 




— 


Rd 




0.22025 




— 


Rs 




0.22025 




— 


Re 




0.17154 




+0.018 


Rb 




0.21585 




-0.440 


AU 

^FB 


-0.03342 


0.06507 


0.12445 


+0.006 


Ad 

^FB 


0.06141 


0.09982 


0.12447 


-0.001 


A^ 
^FB 


0.06150 


0.09983 


0.12908 


-0.001 


A'^ 
^FB 


-0.03350 


0.06513 


0.12451 


+0.103 



Table 6: Summary of Standard Model inputs used in the measurement and their influence 
on the measured At 



FB- 



10.6 QCD correction 

The radiation of gluons from the original quark pair can in principle lower the asymmetry. 
The effect is expected to be small with the present method, where the thrust-defined 
hemisphere charges are calibrated back to the charges of the primary quarks in these 
hemispheres using partly data and partly MC. However, this built-in correction may not be 
sufficient, since the hemispheres defined by the qq axis are not identical to those defined by the 
thrust. Therefore the asymmetry measured by repeating the analysis on the high statistics 
bb MC is compared with the generator level asymmetry in the same selected sample of MC 
events. The generator level asymmetry is found to be higher by a factor 1.0027 ± 0.0009. 

The most accurate calculation available of the full QCD correction in an unbiased sample 
of bb events results in a correction factor 1 + Cqcd 



1.0354 ± 0.0063 [0, g§. The ratio 
of Cqcd to the same quantity calculated with the ALEPH MC is used to scale the residual 
QCD correction. Hereby the residual correction factor becomes 1.0034 ± 0.0019, where the 
error takes into account the statistical uncertainty and the fraction of the theoretical error 
which affects this analysis. 

A summary of all the systematic errors is given in Table |^. 



Flavour purities 

Charge separations 

Detector asymmetries 

QCD correction 



Total systematic error 



±0.00054 
±0.00093 
±0.00016 
±0.00019 



±0.00110 



Table 7: Summary of systematic errors on the Apg measurement. 



11 Cross-checks 



11.1 Simultaneous fit to Rb 

The selection efficiency ffi which yields the flavour purities of the selected event sample has 
been checked by performing the fit on MC and verifying that the input efficiencies are indeed 
returned exactly. It has furthermore been checked that the distribution of 6-tags in the MC 
agrees with that measured in the data after re-weighting the MC by the ratio of ffited to 
simulated flavour efficiencies. 

Finally, the overall bb fraction, Rb, is left free in the fit. In order to get a stable fit, it 
is necessary to tie down some of the previously free parameters, and therefore the small b- 
hemisphere efficiencies for 6-tags less than 0.6 are taken from MC in this case. The result (for 
the three million preselected hadronic Z decays at peak energy) is Rb = 0.2174±0.0008(stat). 



This measured value is consistent with the world average (0.21643±0. 00073 ||2^). The change 
in the asymmetry from letting Ri, float in the fit is —0.0005, in agreement with the sensitivity 
shown in Table 10 and consistent with the uncertainty of the ffi. 



11.2 Lepton tagged events 

For the 6 fit, it has also been checked that it returns the input 6/s when applied to the MC. 
The sample of events with an identified high-p^ lepton offers an additional check on the 
assumptions entering the 6-charge tag. This check is, however, limited by uncertainties in 
the lepton tagging due to c — * /, 6 ^ c ^ / and B^ —>■ B^ transitions as well as lepton mis- 
identification. These uncertainties are reduced to below the 10~^ level by placing a high-p^ 
cut of 1.5 GeV/c on the lepton sample. The distribution of the signed difference between the 
two hemisphere charges in such events is shown in Figure |^. The sign of the lepton provides 
the bb orientation in the event and allows a measurement of 6^ = —0.3395 ± 0.0016, after 
subtracting the 3.8% non-6 background in this sample. This is a factor 0.990 ±0.006 smaller 
than the same quantity measured in the qq MC. The corresponding ratio between the ^^'s 
was found to be 0.995 ± 0.004 in Section p. This then provides a check on the MC prediction 
of kh. 



[l + h 



\MC 



;/- 



^MC/^MC 



0.995 ± 0.007 (stat) 



11.3 Stability checks 



In order to study the stability of the analysis the accepted region in the plane of the two 
6-tags is varied over a wide range. The results for peak data shown in Table |] are consistent 
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Figure 5: The difference between the two hemisphere charges, signed by the lepton charge, 
for events in the E sample containing a high-p^- lepton. 

within statistical errors. 



events selected 


b purity 


raw y4pg 


purity corrected 


charge corrected 


380086 


0.967 


0.0999 


0.1006 


0.1002 ± 0.0031 


491011 


0.938 


0.0987 


0.0996 


0.0994 ± 0.0028 


596618 


0.886 


0.0990 


0.0999 


0.0997 ± 0.0027 


643773 


0.853 


0.0995 


0.1006 


0.1007 ±0.0027 


797132 


0.749 


0.1001 


0.1006 


0.1007 ±0.0027 



Table 8: Results on peak data for various event selections characterised by the sizes of the 
samples and their (predicted) b purities. The event selection shown in bold-face is the chosen 
one, the E-sample. 

The measurements performed for each LEP I year in the data recorded at Z peak energies 
are shown in Table pi. They are seen to be consistent with being constant. 
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year 


^FB 


1991 
1992 
1993 
1994 
1995 


0.0884 ± 0.0125 
0.0992 ± 0.0063 
0.1079 ± 0.0075 
0.0997 ± 0.0040 
0.0966 ± 0.0072 


average 


0.0997 ± 0.0027 



Table 9: Results on peak data for each LEP I year. 

12 Results and conclusions 

The fit to Equation @ is performed in three intervals of collision energy, -^i, with the results 
shown in Figure ^ The asymmetries are then multiplied by the QCD correction factor of 
1.0034 giving, for the sample of events recorded closest to the Z peak, the value: 



^fb(v^ 



91.232 GeV) = 0.1000 ± 0.0027(stat) ± 0.0011(s2/st) . 



The corresponding values averaged over two off-peak energy ranges are 



^fb(v^ 



89.472 GeV) = 0.0436 ± 0.0119(stat) , 
J^^^{^s = 92.950 GeV) = 0.1172 ± 0.0097(stat) . 

The variation of A^-^ with centre-of-mass energy is shown in Figure and compared with 
the Standard Model prediction. The three values are extrapolated to Mz = 91.1874 GeV 
using ZFITTER p^, giving the combined value: 



A'' 



fb 



i^/s = Mz) = 0.0977 ± 0.0025(stat) ± 0.0011 (syst) 

+0.103 (A^B - 0.06513) - O.UO{Rb - 0.21585) 



In order to interpret the result in terms of the Standard Model, a QED and Z — j 



interference correction of 0.0038 is applied according to |23| to give the pole asymmetry, 



.0,6 



y4pg, which is then iterated using BHM 
reached. The result is 



0,fe 



until consistent values for A^'r, A^r and Rh are 



FB5 ^fb 



A 



0,fe 
FB 



0.1009 ± 0.0027 ± 0.0012 



using Ap-g = 0.0623 and Rb = 0.21653. The errors have been obtained from iterating the 
value one standard deviation above the measured ^pB- The corresponding weak mixing 
angle is 



sin^^^ 



0.23193 ±0.00056. 



This is to date the best Ap^ measurement [2,26-34] and provides the single most precise 
measurement of sin^6'^ at LEP. The value is in agreement with the present world average 
value |2^ . 
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Figure 6: The forward-backward bb asymmetry as a function of thrust polar angle at three 
center-of-mass energies. The errors shown are statistical only. 
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Figure 7: The y^ dependence of ApB- '^^^ errors shown are the statistical errors and the 
curve is the ZFITTER prediction, adjusted to reproduce the average fitted asymmetry at 
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